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The interaction potential of a two-dimensional system of excitons with spatially separated 
electron- hole layers is considered in the strong magnetic field limit. The excitons are assumed 
to have free dynamics in the x-y plane, while being constrained or 'polarized' in the z direction. 
The model simulates semiconductor double layer systems under strong magnetic field normal to the 
layers. The residual interaction between excitons exhibits interesting features, arising from the cou- 
pling of the center-of-mass and internal degrees of freedom of the exciton in the magnetic field. This 
coupling induces a dynamical dipole moment proportional to the center-of-mass magnetic moment 
of the exciton. We show the explicit dependence of the inter-exciton potential matrix elements, and 
discuss the underlying physics. The unusual features of the interaction potential would be reflected 
in the collective response and non-equilibrium properties of such system. 
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I. INTRODUCTION 

The availability of new materials and interfaces in semiconductors has allowed the exploration of novel electronic 
systems with fascinating physical behavior. Of particular relevance to the model studied in this paper are the structures 
achieved by clever use of multi-layer geometries, yielding double quantum wells, and heterojunction interfaces of type 
II. In those systems, either by the application of external electric fields or by the intrinsic structure potentials, it is 
possible to achieve separation of electrons and holes into distinct parallel layers, while controlling the in-plane carrier 
densities. 

This situation, of spatially separated electron and hole layers has attracted the attention of several groups, both in. 
theory and experiment. The early proposals of Kogan and Tavger JH as well as Lozovik and Yudson,B and Shcvchenko,cl 
were focused on the possible correlations in such systems due to the electron-hole interactions across the layers. .More 
recently, other authors have theoretically explored different features of these systems, from possible vortices,u and 
dark excitonic states due to hidden symmetries,cl to the various non-trivial thermodynamic phases of these systems.tl 
On the experimental side, there has been substantial activity as well. The experiments of Fukuzawa and coworkers 
gave tantalizing evidence for the—anticipated Bose condensation of 'spatially indirect' excitons in double quantum 
wells under strong electric fieldsE Although later work has shown that the interpretation of those results was not 
reliablejs given the characteristics of the samples used, the concepts of achieving Bose condensation of excitons in 
quantum wells is sound, experimentally feasible, and currently being pursued in new geometries and systemsB 

Controlled electron- hole separation in different layers/planes has also been achieved using heterojunctions of type 
II, such as those formed between InAs and AlSb (or GaSb). In these structures, the band alignments are such that 
electrons and holes are spatially separated in equilibrium, as the bottom of the conduction band on one side of the 
heterojunction lies lower than the top of the valence band on the other. Butov et a/.L2lhave reported photoluminescence 
experiments in AlAs/GaAs heterojunctions, and their results suggest the appearance of a Bose condensate in at least 
the high magnetic field regime. Kono et al. have also reported interesting spectroscopic data suggesting an infrarecL 
active state in the InAs/Al^Gai-^Sb system with unusual properties, reminiscent also of those of a condensate. Ell 
Although the carriers in these latter systems are not introduced optically (as in the experiments above), the close 
proximity of carriers across layers, while remaining at relatively low densities, may yield exciton-like bound states of 
electrons and holes. 

Depending on the details of each system, one can identify suitable conditions under which the electron-hole layers 
would be well described as a collection of polarized exciton-like dipoles .□ These conditions require the in-plane separa- 
tion of charge carriers to be much larger than that across the layers (so that the electron-electron or hole-hole distances 
in each plane are larger than those between electron and hole planes). This in turn yields a system of excitonic dipoles 
predominantly polarized along the normal to the interface. We present here a study of the interactions between the 
resulting exciton states, taking into account both the presence of an intense magnetic field, and the internal structure 
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of the electron-hole pair. The presence of the magnetic field introduces a dynamical coupling between the center of 
mass of the exciton and its relative coordinate, so that the exciton-exciton scattering is a much more complicated 
event than that occurring between point charges. We discuss here these interactions, the potential characteristics and 
different scattering events possible. 

We should mention that perhaps the closest analog of these polarized interacting excitons is that provided by 
polar molecules, such as CO or HF. A large number of theoretical and experimental studies of the scattering events 
between such molecules exists in the literature. 113 Although such systems have permanent dipole moments and live in 
three dimensions, the most different aspect to the excitons here is the non-trivial coupling of their internal degrees of 
freedom with the center-of-mass magnetic momentum. This feature adds a very interesting and subtle complication 
to the excitonic systems studied here. 

In fact, we will show that unlike more compact composite objects, or in the polar-molecule analog, the scattering 
events here can strongly affect the internal state of the participating excitons. In fact, as the in-plane dipole moment 
of the exciton is proportional to its center-of-mass magnetic momentum, the scattering will in general re-orient the 
dipole in a well-defined way which depends on the momenta of the participating excitons. The event may also cause 
transitions to excited internal states of the exciton, just as in the polar molecule analog, although those may be 
suppressed here by the strong magnetic field. This article studies the details of such momentum exchanges and 
effective interactions, and provides explicit expressions for the lowest matrix elements. Apart from describing an 
unusual and interesting situation, these interactions would play a vital role in a description of the collective modes of 
this interacting boson gas. 

It is also interesting to note the similarity of these exciton dipoles with those believed to exist in the quantum 
Hall regime at half-filling of Landau level. The composite fermions there develop a dipole moment proportional to 
the momentum, in a similar way to the excitons we describe. Although the underlying physics is quite different, the 
scattering events of the effective quasi-particles are possibly quite similar J13 Perhaps some of the intuition developed 
in our study of excitons would be of some use in better understanding composite fermions in that regime. 

In what follows, the specifics of the model are described in section O, including a description of the role of magnetic 
field in coupling the various degrees of freedom. Section [II describes then the potential matrix elements for the 
two-exciton scattering events. Section ^ illustrates the resulting scattering potential by considering a few special 
events. Section closes the work with discussion and conclusions. 



II. MODEL 



A. Exciton wavefunctions 



The system of interest can be characterized as a gas of electric dipoles, which are free to move on the x-y plane and 
are effectively polarized either by the application of an electric field in the z direction, or by the built-in heterojunction 
potentials of a type II system, as described in the introduction. For concreteness and simplicity, we shall consider a 
model where the electron and hole layers are separated by a set distance d, and assume that their z-axis dynamics is 
strongly confined. Consequently, the effective layer width of each of the layers is assumed to be so narrow that the 
carriers have only two-dimensional dynamics. This assumption, reinforced by the electron-hole interaction, implies 
that the wave function spread for both electrons and holes in the ^-direction is negligible, and that the other 'transverse' 
states are so high in energy as to be inaccessible for typical situations. We assume further that there is no electron 
tunneling into the hole layer and vice versa. This is in fact the situation for type II heterojunctions due to the built-in 
potentials, and also for indirect excitons in double quantum well systems under strong electric fields. (In the latter, 
however, ihe long-lived indirect excitons co-exist initially with the short-lived direct-excitons created during optical 
pumping.Q) These physical considerations can be suitably represented by constraining the motion of electrons and 
holes to regions z e < 0, and d < Zh, respectively. This approximation .neglects small wave function penetration in 
realistic systems, but given typical parameters, the penetration is smallo 

Although electron and hole cannot overlap in this simplified model of the interface, they still interact via their 
Coulomb attraction, and are able to form a system of spatially separated but bound (if weakly) excitons. As mentioned 
above, this picture of nearly-isolated and well-formed excitons should be an appropriate description whenever the 
inter-particle distance in the plane is much larger (low density) than the electron-hole separation across the interface. 

The wave function for each electron-hole pair in the system may then be written as, 

ip(r e ,z e ;r h ,Zh) = y(r e ,r h )5(z e )S(zh - d), (I) 

where r e and Yh are two- dimensional vectors on the x-y plane for the electron and hole, respectively. This factorization 
makes the problem effectively two-dimensional. Allowing for motion in the z-direction does not alter qualitatively 
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the two-dimensionality, but would require the inclusion of a form factor to account for the finite extension of the 
wavefunction in that direction. This change would only affect the results in a quantitative way, and can be introduced 
straightforwardly, as done before in similar situations. 

The magnetic field is applied perpendicular to the layers, H = Hz, and it is assumed to be sufficiently strong that 
the relation 



r H <^a e ,a h (2) 

holds, where a e< h — n 2 n/m e .he 2 are the effective Bohr radii for the electron and hole, m e ,/i are the effective masses at 
H = 0, k is the average background dielectric constant, ejs the charge quantum, and rjj = y/hc/eH is the magnetic 
length. As discussed first by Gorkov and Dzyaloshinskii,E_3 condition (||) allows one to apply perturbation theory in 
this rather complex problem. Motion of a 2D neutral electron-hole pair in a transverse magnetic field H = (0, 0, H) is 
described by a Hamiltonian, which for non-degenerate and isotropic bands in the effective mass approximation reads 

H = ilHW e + + J-HW, - lA h f - . (3) 

Here, n = (/ci + k,2)/2 is the average background dielectric constant across the heterostructure. Since Ki ~ K2 in 
typical systems, possible image charge effects are small and neglected here. 

The dynamics of the single exciton system is characterized by a conserved quantity associated with the operator 
for magnetic momentum of the center of mass, P = — ihVn — (e/c)A(r)0 Here, 

R = (m e r e + m h T h )/M (4) 

is the center of mass coordinate, while 

r = r e - r h (5) 

is the relative coordinate of the electron-hole (e-h) pair, M = m e + m^, and A(r) = |H x r in the Landau gauge. 
The wavefunctions $„ m p, describing the state of an e-h pair in the field H can be written as,E3 



*™p(r e ,r ft ) = exp j - (P + — H x r J • R 

exp Q^7 p ' r ) *„m(r-rp), (6) 

where 7 = (m/, — m e )/M, P is the "center-of-mass" or "magnetic momentum" of the exciton associated with the 
operator P, 

r? 



rp = f«P, (7) 



and the wave function $ nm is identical to the wavefunction of a charge e in a field H 







n 



I 1 / 2 „- 



e im<p 



X 



2|m|+i( n+ I m |)!7r 



*^nm (-0 

/ \ l m l 

where are Laguerre polynomials, and p =| r |. The wavefunctions \&nrrfPi(r e , ^h) describe then, in the limit of high 
magnetic field, the dynamics of magnetoexcitons with dispersion relationJLj 

with 

£ nm = huj H ( n + i J (| m | -7m + 1) , (10) 
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where the cyclotron frequency ujh — eH/fxc is defined in terms of the reduced mass of the exciton /z, and 



E nm {P) = ($„ m | ; ; r | $„ m ). (11) 

K | r + rp | 

Equation (g) is the energy of the $ mra p state to first order in the Coulomb interaction, and is a good approximation as 
long as (||) holds. The states constructed in this fashion can be viewed as an exciton that has center-of-mass motion 
P in the x-y plane, and with 'internal structure' given by the state <i> nm (r — rp). 

It is important to emphasize that the functions $ nm in (^3j) are centered at rp, given by Eq. (7), so that the actual 
in-plane separation between electron and hole is (r) = rp, proportional to its magnetic momentum. Notice moreover 
that the average radius vector between electron and hole is orthogonal to P, for all n an m, and it vanishes for P = 0. 
This peculiar dependence of the dipole moment on the magnetic momentum can be intuitively understood as the 
result of the Lorentz force tending to separate the charges in each pair. Notice that the in-plane polarization reduces 
their Coulomb interaction, and makes the exciton more susceptible to ionization by system imperfections .113 This 
P-dependence will have also important consequences for the scattering matrix elements, as we will see later: even 
an elastic re-orientation of P results in a re-alignment of the dipole moment, which in turns changes the interaction 
with other dipoles. Nevertheless, notice that because this system has spatially separated electron-hole layers, the total 
dipole moment vector of the exciton has a constant component along the z-axis, and this is the dominant component 
in most cases. 



B. Inter-exciton residual potential 



As already mentioned, the charge separation imposed by the layered geometry produces an effective polarization, 
nearly perpendicular to the interface for small P values, and resulting in a non-zero dipole moment for all the excitons 
described here. This fact gives rise to an overall repulsive interaction between all excitons in the system. It will be 
this "residual potential" that provides for a collective response, as we describe in the next section. 

To the lowest order in a multipole expansion, the residual interaction potential between two excitons located at x 
and x', respectively, can be written as the interaction between two dipoles, 

v _ P P' 3[p-(x-x')][p'-(x-x')] 

\ I 13 I MB ' V / 

K | X — X' |° K I X — X' |° 

The dipole moments p are generated by the non-zero expectation value of the relative coordinate. Correspondingly, 
p = e(r + zd), where r is the in-plane relative coordinate, and d is the z-axis separation. As this expression depends 
on both the relative and the center-of-mass coordinates (x and x') it needs to be evaluated for each exciton state 
wavefunction (see next section). 

Notice that this dipolar approximation should be valid as long as the exciton separation is larger than any of the 
characteristic size-length scales of the excitons themselves, 

|x-x'|>|(r)|,d,a M , (13) 

where (r) is the in-plane exciton mean radius (= rp), a M = Ti 2 nj fie 2 is the exciton's Bohr radius, and d is the z-axis 
e-h separation. For closer inter-exciton separations, one should in principle include higher multipoles in the interaction 
between excitons, accounting for the constituent electrons and holes. The dipolar approximation would break down 
as the in-plane carrier density increases, violating the condition (|l3|). Correspondingly, this condition would require 
|x — x'| > a M 100 A in typical materials/systems. The in-plane densities would need to satisfy n.^; l/^a 2 « 10 12 
enr 2 , quite a reasonable request (given the typical experimental densities of 10 10 to 10 11 cm _2 )J3illl 



III. INTERACTION MATRIX ELEMENTS 



As discussed earlier, the main motivation for considering this problem is to model the dielectric response function 
for a 2D Bose gas of dipole-like polarizable bosons in a strong magnetic-fcld. In this context, the dielectric function 
in the self-consistent field mean-field approximation can be written asj 



tic. iei< 

Va'PiaP'HpP' ( w ) 



(14) 
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where Upp' is the polarization matrix, u> is the frequency of the perturbing potential, and the inter-exciton interaction 
potential matrix elements in the excitonic wave function basis are given by 



r al mUT')v(T,T')^ a (r)^,(r')dTdT'. 



(15) 



Here, T and V refer to the exciton degrees of freedom, with T = {r,R, z e ,Zh}, or {r e ,rh, z e , Zf t }, and the a and (3 
indices denote the {nmP} set of excitonic quantum labels. Using the states described above, we have for the first 
term of the potential (12), 



while for the second term, 



v a'(3;uf3' 



Pa'a 



(16) 



y(2) 



a'f3;a(3' 



(q-Pa'a)(q- P/3/3' )Xl{l) 
"(Pa'a ' P/3/3' )X2(<?) <5(q-q') 



(17) 



The delta functions in these equations ensure overall magnetic momentum conservation in the scattering between two 
excitons, i.e., 



P' + K = P + K' , 



(18) 



with Tiq = P — P', and fiq' = K — K', as one expects for translational invariant systems. The labels for incoming 
and outgoing momenta and other exciton quantum numbers are shown in Fig. [l], where q represents the in-plane 
momentum transfer due to the scattering event between excitons. In this notation, the scattering process is fully 
described by the change of the remaining internal state labels, a' — > a, and b — > b' , as indicated in the figure. 
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FIG. 1. A sketch of a scattering event between exciton dipoles. One of the exciton-particles makes the transition 
{P',a'} — > {P,a}, while the other one changes {K, b} —* {K',fe'}. Notice that in-plane dipole moments are perpendicular 
to each momentum, and they change upon scattering (as indicated by blank and shaded arrows). Dashed line indicates the 
momentum exchange due to the inter-exciton potential. 



In these expressions, the functions 4>(q), Xi(q)i an d X2(<z)j are Fourier transforms of the dipolar interaction depen- 
dence on the center-of-mass coordinates, and are explicitly written in Appendix A. The in-plane "dipole moment" 
matrix elements are given formally by, 

2e d 

Pa'a = —-^-M a > a , (19) 

i"f oq 

where Tiq = P — P' as above, and the non-local "overlap" matrix elements are given by 

= e*** 1 *" J e^ r S;, m , (r) $ nm (r - u) d 2 r , (20) 

with similar expressions for ppp> and Mpp> . In the last equation, we have used u = rp — rp/, and the gradient in Eq. 
( pj| ) refers to the explicit g-dependence shown in (pfj|), different from u. 

In the following sections, we describe typical features of the overlap and dipole matrices, M aa i and p aa ' ■ 



A. Dipole matrix elements 



In order to better understand the nature of the residual interaction matrices, we evaluate some of the lowest 
elements. In what follows, and for notational convenience, we use a = {nmP} — {aP}, with a standing for the 
indices of the <& n m states, so that one can write for example, 

Pa'a = Pa P .aP = Pa>a(P ~ P', P + P') = Pa'a(q, P') , (21) 

where the explicit dependence on the sum and difference of participating momenta is indicated. The last equality 
uses the overall conservation of magnetic momentum provided by the delta functions in ( |l6| ) and Jl7|). 
The simplest dipole moment matrix element (for a — {00} = a') can be written as (see Appendix B), 

Poop'.oop - Poo,oo(q, P') = | (r 2 P'+ fiq + t^q) e***'"--^ e -d+7 2 ) 9 2 ^/8 ; (22) 

where the in-plane momentum exchange q, and the incoming magnetic momentum P' are used to specify the dipole 
and non-local overlap matrix elements. 
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From this expression, the (non- interacting) long- wavelength limit q = yields, 



Poo,oo(,q 



0,P') 



2 
erp 



r 2 p- = [z x P'] 



(23) 



This represents what one could call the 'proper' dipole moment of the exciton in state nm = 00, and with magnetic 
momentum P', since the expectation value of the relative coordinate is rp< . In fact, it is possible to show (see Appendix 
B) that all the diagonal dipole matrix elements in the limit q — > yield p a a(q = 0, P) = erp, since in fact all <& nm 
states have the same dipole moment (in this high field limit), as we discussed following Eq. (|l l|) . We emphasize that 
large momentum values correspond to larger exciton size and lower binding energy, as the exciton is increasingly 
polarizedJlj One expects that such high-P states would be easily affected (even disintegrated) by perturbations in the 
system, such as impurities and surface inhomogeneities. 

It is interesting to note the role that 7 = (m/j — m e )/M plays in Eq. (p2|), providing an imaginary part (or phase) 
proportional to q to the dipole matrix element. Notice further that for any 7 values, a non-vanishing momentum 
exchange q depresses exponentially the dipole matrix element, with a characteristic length ss ru- Since non- vanishing 
q corresponds to the momentum/energy transfer from one exciton to the other, high momentum transfer processes 
will appear to be strongly suppressed by this potential. Let us discuss these features in the next section. 

B. Potential matrix elements 

The simplest elements of the potential are those diagonal in the {a, a'} indices. For two excitons with incoming 
momenta P' and K which exchange momentum q, the potential matrix element is given by, 



Vr 



00,00;00,00 



(q,P',K) = U(<z) 



-(r 2 p'+7iq + hr 2 H q) ■ (r 2 K-hq - ijr 2 H q) + e 2 d 2 



V, 



00,00;00,00 



(q,P',K)| x 



M 00 ,oo(-q, P')M)o,oo(q,K) 



where, 



and, 



M) ,oo(q, P') = exp ( i-q • r 2P < +?iq - (7 2 + 1) 



r 2 ci 2 



(24) 



(25) 



2 

Voo,oo : oo.oo(q, P > K ) = -[{ °l ( r 2P'+fiq + i7r 2 H q 

-(■ 



q- I r 2 K-hq - ijr H q 



r2P'+?iq + tjr H q ■ I r 2K -?»q - H r H 



r q)x2(<z)} ■ 



(26) 



Notice that these expressions contain a contribution from the constant ^-component of the dipole moment, ed, as well 
as from the in-plane components. 

From these equations, and considering the q — > limit of the potential (see Appendix A), we may write 



Vr 



00,00;00,00 



2 

(q -> 0, P', K) = { (l - 6 ^' 2 r(3/4)) r P , • r K + d 2 } 



(27) 



This result is expressed in terms otthe proper dipole moment of each exciton, proportional to rp< and tk, as intuitively 
expected by Lerner and Lozovik.ll3 It is clear that the sign of the first term in the interaction (|27j) depends on the 
relative orientations of P' and K, and the resulting dipole moments. The total interaction between excitons will 
be more/less repulsive for anti/parallel P' and K, as the contribution to the z-axis moment is modulated by the 
in-plane component. For small and moderate magnetic moment values, typical of the excitons in this system at 
low temperatures, the repulsive interaction is however only weakly modulated, since d > r 2 p = r 2 H P, but it is still 
dependent on the relative orientation of the proper dipoles. 
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IV. SCATTERING EVENTS 



As we have mentioned, the potential matrix elements above allow the description of the collective excitations of 
the weakly repulsive gas of polarized excitons, in a manner similar to that treated in Ref. [|l^,|2^]. Moreover, these 
potential expressions can also be used in a quantitative description of the kinematics of scattering events, as those 
needed in a treatment of the distribution function via the Boltzmaim equation to evaluate drag,Ell or the evaluation 
of the Bose condensate properties in this dipole-interacting system.E3 As an illustration of their use, we describe in 
this section how the potential matrix elements calculated above can provide rates and cross sections for different 
inter-exciton scattering events. For simplicity, we deal here with 'elastic' collisions, when there is no change of the 
internal state under scattering. More complex events are in principle also allowed, although 'inelastic' processes are 
suppressed by the strong field. 

If one considers scattering events in which the internal state of the excitons is left unchanged, one is then faced 
with a purely kinematically elastic collision. The description of this elastic collision, like any problem of two bodies, 
is simplified by changing to a system of coordinates in which the center of mass of the two particles is at rest. The 
scattering angle in the center of mass reference frame is denoted by 9, and it is related to the angles 9\ and 02 giving 
the scattering angles of the two particles in the laboratory system-of coordinates. In the case in which the second 
particle was at rest before the collision, for example, one can writetll 

m 2 sin 6 1 . n . 

tan0i = f -, 2 = -tt-0, 

mi + m 2 cos# 2 

where mi, m 2 are the masses of the scattering objects. In our case, the masses of the two 'particles' (excitons) are 
the same (mi = m 2 = M^ xc , which depends on the internal state a — {nm}), and we have simply, 

*i = ^, e 2 = l -^-e)- 

so that the particles diverge at right angles in the laboratory frame. 

The scattering cross-section can be calculated using the Born approximation, since the residual potential (p^), may 
be considered a weak perturbation. Notice, furthermore, that the residual potential depends only on the distance 
between excitons x — x', so that the scattering field is central. Now, in the center of mass frame of reference, we can 
write 

P' CM = n K CM = -n , (28) 

where n = (P' — K)/2, is the relative magnetic momentum between excitons. Thus, the interaction potential matrix 
element for this event, where the internal state is assumed to be a = {00} before and after the collision, can be written 
as 

^oo,oo ; oo,oo(qi n ) = \ j 4 <M<?)d 2 - X2(q) ~ <j>(q) (r 2 „+Kq + ijr 2 H q) - (q • r 2 „ + i-yr 2 H q 2 ^ Xi(<?) j x 



r 2 a 2 

-(7 2 + l)^ 



(29) 



where fiq = P P' = K — K'. Other internal states are given by a different detailed expression, but identical 
kinematics (Appendix B). 

In the case under consideration, Eq. ( p9| ) describes the matrix element for a 'transition' (scattering event) from a 
state with momentum n to the state with momentum n' = (P K')/2 = n + fiq, which we could then denote as 
U nn >. Correspondingly, the scattering rate can be calculated from the golden rule, 

dW fi = {2Tr/h)\U nn ,\ 2 5(E' n ~ E n ) , (30) 

where the final and initial energies of the exciton of interest are, 

E' n -E n = (n' 2 -n 2 )/2M^ c , (31) 

where Mqq c is the effective excitonic mass of state {00}, including the electron-hole interaction, for |n|, |n'| -C h/ru .0 
The presence of the delta function in (|3(]) indicates that the scattering event is kinematically elastic. 



8 



V. CONCLUSIONS 



We have presented explicit expressions for the residual interaction between polarized excitons in strong magnetic 
fields. This potential would have important consequences on the description of individual scattering events, as dis- 
cussed, and the collective excitations of the system. Unfortunately, it is not clear how one would perform experiments 
to directly measure the many different scattering processes possible. Although we are hopeful that some experiments 
might be designed in the future to analyze these processes, we believe that the more direct probe would be study of 
the various collective modes in this interesting and unusual system. Since the density-fluctuation modes are now able 
to include rather complex internal excitations, the resulting modes may indeed be quite unusual and complicated. 
One would use the potential derived here in an approach similar to the case with no fieldJi3'E3 and results will be 
presented elsewhere. 

We trust that these expressions would also be useful in the description of other kinematic and thermodynamic 
properties of the system. 
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APPENDIX A: POTENTIAL FUNCTIONS 



Because of the spurious short-range divergence introduced by the |x — x'| 3 dependence on the dinolar potential, 
we use the physical cutoff parameter provided by the finite z-direction polarization of the exciton d.cB Therefore, by 
introducing the regularization factor as per |x — x'|~ 3 — > |x — x'|~ 3 (l — e~l x ~ x I / d ), we preserve the long-range 
dipolar interaction, while allowing for the short-range Coulomb-like repulsion. Correspondingly, we may write, with 
S = x — x', 



4>{q) = ~j e^a-e-^Vs 



2ttk 4 



q 2 d 2 \ / 4 \ (q 2 d 2 



(Al) 



where /„ is the Bessel function of imaginary argument. Notice, incidentally, that 



4>(q - 0) = ^-j-p , (A2) 

and is therefore ill-defined for d = 0. This function 4>(q) accounts for the first term in the interaction in Eq. ( |l7|) or 
(§)■ . 

Similarly, the second dipolar term can be written as 



/(2) = i /• e i q -s (p-s)y-s) (i _ e _ sVd , )d2s 

K J D 



_ S" 



where for qd 3> 1 



and 



(p-q)(p'-q)Xi(?)-X2(g)p-p', (A3) 



3^ 27r(-256) fc (4fc + l)(4fc + 3)r 2 (fc + l)r 2 (fc + |) 

Xl(q) r 2 (i)r 2 (f)(£: + i)w 4fe +v fe + 5 1 ' 



3 ^ 27r(-256) fc (4fc + i)r 2 (fc + i)r 2 (fc + 1) 
K t^ r2(i)r 2 (|)(fc + i)!^+v ,!+3 



On the other hand, for qd <C 1, Eq. (A3) behaves as 



f(a?= 3^r(3/4) p 



(A6) 



APPENDIX B: DIPOLE MATRIX ELEMENTS 



Notice that the dipole moment matrix elements obey the symmetry relation 

# 

PnmP,ii'm'P' — Pn'm'P' ,nmP • 

Some special cases of overlap and dipole moment matrix elements follow. For the diagonal elements, 



(Bl) 



Moop', oop = exp 



where u = rp — rp< , and 



ilq-lrp+v)-' 



,2 



2 2 

Poop<,oop - I (rp + r P , + ir^yq) e ^<i-(«-+rp')-(7 3 +l)-V- 
= 2 (rp + r P / + ir%jq) Mqop'.oop 



The si 



dmplest/lowest off-diagonal elements are 



and 



Similarly, 



2 r 2 2 u 2 \ 

!- T ux q-z Mqop'.oop 



V 8 8r 2 H ' 4 



P01P',01P — 



1 - r ^ 2 ( ^- + ~ + Poop'.oop 4- »(7 + 1)^^Moop',oqp 



Mn-iP' 0— IP = 1 — rr~ — — 

l « Srj r 4 

\ 2 2 (l 2 , 7 . 1 

1-r * 9 U + I + 8 



7 2 r| f g 2 w 2 7 A 
8 8^ 2__ 4 UXq ' Z J M ° 0P '> 00P 



Po-ip'.o-ip — 



1 /. 7 r H qe- 1 ^ ue-^ 

J«noP' oi p = — 1= * 

oop ,oip ^- ^ 2 2r ff 



er^q ~ 

Poop',oop + H7 — 1) — 2 — ^oop',oop ; 

AfoOP',OOP j 



Poop',oip = i ~2~^ TH t q ' - *y)]PooP',oop + T/l"^ - *y) M ooP',oop ; 



1 ( .ir H qe lv " ue^ 



A/oOP',OOP 



Moop',o-ip = ^— 

'Y — 1 &T*H 
P00P',0-1P = r H [q ' ( x + »y)]P00P',00P + T/l"^ + *y) M 00P',00P 



AfoiP',0-lP — 

Poip',o-ip = — {j 2 + 1 



7 2 r 2 ^ g 2 e i2¥, '' u 2 e i2v ^\ 

" 1 I •'"OOP', OOP : 



8r 2 



H 



.r 2 H q 2 e l2v i .er 2 H qe liPq ,„ ... - 

~ g Poop',oop+« — -y= — (x 4- iy)M 00 p',oop ■ 



(B2) 



(B3) 



(B4) 



(B5) 



(B6) 



(B7) 



(B8) 



(B9) 
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